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Abstract 

We calculate 0{as) corrections to inclusive and isolated double prompt pho- 
ton production, both for the unpolarized case, and for longitudinal polar- 
ization of the incoming hadrons. The calculation is performed using purely 
analytical techniques for the inclusive case, and a combination of analytical 
and Monte Carlo techniques to perform the phase space integration in the iso- 
lated case. A brief phenomenological study is made of the process pp — > 77^ 
at CMS energies appropriate for the RHIC heavy ion collider. 
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I. INTRODUCTION 



Among the various and still unexplained aspects of QCD, the study of the spin structure 
of the nucleon has gathered considerable work in recent years , suggesting that a simple 
parton model interpretation of the phenomenon is far too simplified to completely describe 
the distribution of spin inside the nucleons. The original EMC data suggested that quark spin 
accounts for only a small fraction of the nucleon spin. Since then the possibility that gluons 
may carry a significant fraction of the nucleon's spin has been advanced and investigated 
theoretically (see ^ and Refs. therein). But other alternative explanations, such as the 
possibility that the proton may have a large negative strange sea polarization have also 
been suggested, along with various scenarios which use a combination of the two extremes. 

Recently, new experimental results on the polarized proton structure function have 
become available [Q. Most of the information collected so far from the phenomenological 
side, which however still leaves an appreciable disagreement with naive parton model expec- 
tations, is based on studies of the hadronic tensor in the kinematic region of Deep Inelastic 
Scattering (DIS), through its structure functions gf^x^Q"^) (x is the Bjorken variable), pri- 
marily at large and finite x. This discrepancy, referred to by some as "the spin crisis", 
indicates that sum rule predictions based on the constituent quark model picture of the 
nucleon, need to be amended, since it is only in the DIS limit that the structure functions - 
which parametrize the hadronic tensor - scale to quark distributions which are functions of 
X only. 

At finite Q^, DIS sum rules are violated by quark-gluon interactions, usually parameter- 
ized by higher twist effects in the Operator Product Expansion. The contribution coming 
from the axial anomaly (see also |^ and Refs. therein), in particular, has been shown 
to induce a cancellation between quark and gluon contributions to the first moment of the 
polarized structure function, gi{x). 

Although the anomalous contribution is a radiative effect {0{as)), renormalization group 
arguments and the log((5^) increase of the nucleon spin suggest that such contribution (which 
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is as((5^)AG) remains finite even in the — > oo limit in which the O.P.E. is apphcable. The 
"spin crisis" is therefore solved by the observation that the EMC result is not a measure of the 
quark spin, according to the naive parton model picture. It only measures the combination 

9 V 2n J 

between the quark AS and the gluon contributions to the total spin. This leaves open the 
question as to the actual size of AG, i.e. whether it is large or small. This question can 
only be settled conclusively by experiments which are capable of directly probing AG. 

Given the universality of the structure functions within QCD, it is expected, in a few 
years from now, that the information gained from DIS polarized scattering will be supported 
and supplemented by new results from polarized proton-proton collisions at the BNL heavy 
ion collider RHIC in forthcoming experiments. It is therefore important and interesting to 
see how various processes at large pt are affected -in the polarized case- in their behaviour 
when next-to-leading (NLO) order corrections are included. It is also interesting to see how 
the asymmetry between the unpolarized and the polarized cross sections behaves when a 
reduced scale dependence is present, due to firmer NLO perturbative calculations. From the 
theoretical side, further motivation to proceed with these studies comes from the fact that 
the splitting functions for a complete NLO Q^-evolution of the spin structure function gi 
have been recently obtained 0, and new parametrizations in NLO will soon be available for 
the first time. Many processes have been suggested for study at RHIC 0, but only a few 
so far have been calculated at NLO 0,^. 

In this work we present a next-to-leading order study of the process pp — > 77X, which 
has been calculated in the past [|10| for unpolarized photons at NLO. We repeat both these 



calculations and extend them to the polarized case, which to our knowledge has only been 
studied in LO |Tl|, with inclusion of the higher order (0(a^)) process gg 77. There it 



was found that the cross section should be large enough to measure, and should have some 
sensitivity to the polarized parton distributions but that the asymmetries are not very large 
except in regions where sea quark spin dependent processes dominate. This meant that 



the process would not be useful for extracting information on AG. But they also made 
the point that a full NLO study including all relevant contributions would be necessary to 
confirm these conclusions, since cancellations between the various processes may occur, and 
the effect of the corrections were completely unknown. For example, they could not include 
contributions for the higher order process qg 775', which are known to be quite significant, 
and can even be larger than the qq scattering process in proton-proton collision processes, 
as opposed to proton-anti-proton collision processes. Also, in addition to the possibility 
that double prompt photon production may yield important information on the polarized 
proton structure function, particularly the possibility that a large fraction of proton spin is 
carried by gluons, or a large negatively polarized strange sea, the process is also important 
because it is a background to Higgs decay. A full NLO analysis using more modern polarized 
parton distribution is now necessary to either support or disagree with the conclusions of 



Ref. [jTT|. In this study we are not yet able to use parton distributions evolved in NLO 
as these are not yet ready, but in a future study we hope to include them [|12|. We can 
nevertheless extend that of Ref. |]TT|, by including all 0(as) higher order corrections to the 
hard scattering process and correctly taking isolation effects into account. 

The organization of the paper is as follows. In section II we describe the methods involved 
in the calculation and the regularization prescriptions adopted for the inclusive case. The 
main contribution of this paper is the calculation of the polarized cross section, but since 
the unpolarized cross section has been presented before, we discuss this first in each case 
and then extend the arguments to the polarized case. In section III we outline the combined 
analytic/Monte Carlo method of calculating the isolated cross section. In section IV we 
present and discuss some numerical results, and in section V we draw some conclusions. 

II. INCLUSIVE DOUBLE PHOTON PRODUCTION 

The inclusive cross section we calculate is 

da 



dkj^^dyidz 
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This cross section is differential in the transverse momentum, kj^^, and rapidity of one of the 
photons (referred to as the trigger photon) plus the variable z defined by 

^ = — urTir- (2.1) 

z thus contains information about the transverse momentum of the second photon, but not 
its rapidity The cross section is thus limited in that experimental cuts on the rapidity of the 
second photon and isolation cuts on either photon cannot be implemented. It is nevertheless 
still worth calculating, since it provides a check on the more versatile calculation using 
combined analytic and Monte Carlo techniques which we describe in section III. 

A. LO Contributions 

In LO 0{al^), the contributions to double prompt photon production are from qq an- 
nihilation to two photons, and the various single and double fragmentation contributions, 
some of which are shown in Fig.l. Prom Eq.(2.1) we can see that the qq annihilation process 
in LO is proportional to 5(1 — z), as the k'j^s of the two photons must balance. The hard 
subprocess cross section is 

da , _ , 27rQ;^„ ^l-2v + 2v'^ 



{qq - 77) = -^e^ ^ 5{1 - z)5{l - w) (2.2) 



dvdwdz 3s ^ v{l — v) 

where 



t 

V = 1 + - 
s 



ctem is the electromagnetic coupling constant and denotes the quark charge. The usual 
Mandelstam invariants are defined in terms of the momenta of the two incoming hadrons 
Pa and Pb, the momentum fractions of the initial partons Xi and ^2, and the momentum 
of the trigger photon ki via. 
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S = {xiPa + X2PbY = X1X2S 

t = {xiPa - kif 

u = {x^Pb - hf (2.4) 

where is the center-of-mass energy in the hadronic system. 

In this study we do not include contributions where both photons are produced by 
fragmentation off a final state parton as these are expected to give a very small contribu- 
tion, especially after isolation effects are considered. We do include process where one of 
the photons is produced directly and the other by fragmentation, referred to as the single 
fragmentation contributions (Fig. lb). These processes contribute to O(ag^), although the 
hard subprocess is 0{asOtem)i since the fragmentation function of parton i into a photon, 
D^/i{z,Q'^) is 0(as/aem)- In a fully consistent NLO calculation we should also include 
the 0{as) corrections to these fragmentation processes, but we are not yet in a position to 
provide these. 

The full expression for the physical qq annihilation cross section in LO is 
da 27ikTi dv ,.2N.i?/ ,^2^ da 



, , c / 1 /,(^i,M^)/,(a^2,M^)— — -(gg-^77) + (g^g), (2.5) 

dkjnidyidz irb Jvw 1 — v ^ ^ dvdwdz 

where V and W are defined similarly to v and w in Eq.(2.3) except now in the hadronic 
system, and //^(x, M^) is the parton distribution function for parton of type i in hadron A 
as a function of x, the momentum fraction, and the scale, M^. 

There are two types of single fragmentation contributions to the process. Those where 
the trigger photon 71 is produced by fragmentation and those where 72 is. These contribute 
in different regions of z, the former in the region z > 1 and the latter in the region z < 1. 
In terms of the hard subprocess cross sections for the processes qq — > 'yg and qg 7g which 
can be found, for example in Ref. [0, these contributions are given by the expressions; 

dkridyidz irb Ji-v+vw z'^ Jvw 1 — v dv 

D,^,i{z',Ml)e{z-m-^-z'), (2.6) 



and 



do- , 1 -^r^ /"^ dz' dv ,,9x.R. ^ ^'?^da'^^''^'' 

= 27rkTi— E / — / 1 f^ i^u M')ff{x2, — ^ 

dlvTidyidz irb ~^Ji-v+vw z'^ Jvw 1 — v dv 



D^,/u{z',Ml)e{l-z)5{z-z'), (2.7) 

where i,j and k run over the various quark flavours and the gluon and M|. is the fragmen- 
tation scale. 

In the polarized case the cross sections are given by the same expressions as above, but 
now the hard subprocess cross sections must be replaced by the corresponding polarized ones, 
dAa / dv. The parton distributions must also be replaced by the corresponding polarized ones 
which are defined by 

AA^(x,Q2) = /^+(x,Q2) (2.8) 

where /^_|_(x,(5^) is the distribution of parton of type a with positive (+) or negative (— ) 
helicity in hadron A. Likewise the polarized subprocess cross sections are defined by 

d/S.a da , ^ da , ^ , ^ 

Again, +, — denote the helicities of the incoming partons. Since only the initial state 
is polarized, we do not need to change the photon fragmentation functions. The hard 
subprocess cross sections in LO for these can be found, for example, in Ref. 0. 



B. The Box Diagram 

At collider energies, the box diagram process, gg — > 77 (Fig.lc), although it is ©(a^ag^), 
has been shown to contribute significantly to the cross section for unpolarized double prompt 
photon production in the low kr region due to the large gluon distribution at low-x. It is 
questionable whether or not one should include contributions from the box diagram in an 
0{as) NLO calculation, since it is of higher order, and may introduce numerical cancellations 
with the genuine NLO contributions leading to misleading results. In a proper phenomeno- 
logical study at NLO we would not include this contribution, but in order to make contact 
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with the work of Ref. [12], and to estimate the contribution from this process, we will include 
it. The unpolarized hard subprocess has been used many times before and can be found, 
for example, in Ref. ||13||. For the polarized case the hard subprocess cross section can be 
obtained from the amplitudes according to the prescription given in Ref. |11|]. To calculate 
the contribution form this process, we use a similar expression to Eq.(2.5), with the quark 
distributions replaced by gluon distributions and the hard subprocess cross section replace 
by that for the gg 77 one. 



C. The NLO contributions 

At NLO O(aemtts) there are the virtual gluon corrections to the LO qq annihilation 
process, plus real gluon emission corrections qq — > 77(7 (Fig. (2)). In addition, the 3-body 
process qg 77^ also contributes (Fig(3)). The virtual contributions can be obtained 
directly from Ref. 0, where they were calculated for single prompt photon production in 
both the polarized and unpolarized cases, by simply removing the non-abelian couplings. 
The 3-body matrix elements can be obtained in a similar way. 

A note concerning the matrix elements for the polarized case is in order here. When 
taking the traces to calculate the helicity dependent matrix elements, one needs to project 
onto definite helicity states, labelled h for quarks and A for gluons, for the incoming particles. 
This is achieved for quarks by defining spinors of definite helicity according to 



h) = \r{pi)^{l - /i75) (2.10) 



where p is the momentum. For gluons the general expression is 

1 



£^(^2, A)e*(p2, A) 2 



P a 



(2.11) 



This introduces the 75 matrix and the antisymmetric tensor e^i^po- which are defined in 4- 
dimensions. To work inn ^ 4-dimensions, some consistent scheme must be chosen in which 
to treat these 4-dimensional objects. In Ref. P], the original scheme of t'Hooft and Veltman 
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1^ , and Breitenlohner and Maison [|T5| was chosen. Briefly, this scheme makes a division of 
ra-dimensional Minkowski space into a 4-dimensional and a {n — 4)-dimensional part. Thus 
an arbitrary vector p will have a 4-dimensional part, p, and an (n — 4)-dimensional part, 
p. This means that scalar products of these 'hat' momenta will inevitably be present when 
traces are taken. These terms in the matrix elements must be treated properly when phase 
space integrals are performed. We shall consider this point in detail in appendices B and C. 

In order to obtain a cross section differential in kxi, yi and z, we need to integrate 
3-body phase space over the kinematic variables of the unobserved quark or gluon. The 



method of performing these integrals have been detailed in Ref. |]I0[. We outline the method 
again in Appendix B and extend it to the polarized case where 'hat' momenta must also 
be considered. As in Ref. we restrict ourselves to the cases where the photons are 
in opposite hemispheres, i.e. z > 0. The phase space integrals are performed in 4 — 2e 
dimensions in order to expose soft and collinear singularities as poles in e. 

In the case of qq scattering, we obtain single and double poles in e from the 3-body 
phase space integrals. When the virtual contributions are added to this the double poles 
automatically cancel and the remaining single poles must be absorbed into the parton dis- 
tribution functions. Thus there are three pieces to our calculation for this process which 
must be combined in order to get a finite partonic subprocess cross section. There is the 
virtual contribution which we represent by the function 

^ 1 1 

' S, f , z. 



dvdwdz V ' ' ' ' e 



There is the result of integrating the 3-body matrix element over the phase space of the 
unobserved partons as described in appendices B and C, which we denote by the function 



1 1 



qq y "^J ^) ^) ^5 ^2 ' ^ 



and there is the factorization counter term 



1 da^g _ a. 



sv dvdwdz 271 



1 / 1-v y\ da'i'i^^^ f l-v . , 

H,, , — S, vw, e 5 1 - z) 

1 — vw V 1 — vw J dv V 1 — vw 



1 da'^'^^"''^ 

-Hgg{w, M^)—- {ws, v, e)S{l - z) 

V dv 



(2.12) 



where 



:P,Az) 



Q2 } ^J^3 



(2.13) 



In the MS factorization scheme which we adopt, 1/e = 1/e — 7^; + ln47r, fij{z) = 0, and 
Pij{z) are the well known one-loop splitting functions for parton j into parton i |T^. The 
finite subprocess cross section is now obtained by adding these three parts, 



Kgg{S,V,W,Z,M'^] 



da 
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dvdwdz V ' ' ' ' e 



qq 



1 d^qq 1 
sv dvdwdz V e 



(2.14) 



where is the factorization scale. 

A similar procedure is followed for the qg initiated process, except now there are now 
no virtual contributions and hence no double poles. Also different singularities are encoun- 
tered when integrating the 3-body matrix elements, hence we have a different factorization 
formula; 



1 da, 



qg 



sv dvdwdz 



'271 



V + vw 



H^q ( 1 



V + vw, Mi 



da^'i^^'i 
dv 



V 



V + vw ' 



1 da^^^^'^ 

-Hqgiw, M^)—- {ws, V, e)6{l - z) 

V dv 

+ -H^g (z, Ml) (s, V, t) 5(1 - w)e{l - z) 

v ^ ' dv 



e ) b{zx - z) 



(2.15) 



where 



E^g[x,Ml) 



\P,q{^) 



Ml, 



(2.16) 



Mp is the fragmentation scale and zi = — v + vw). 

Once {A)Kqq and {A)Kqg have been calculated then we can use them to calculate the 
physical cross sections by convoluting them with the appropriate parton distribution func- 
tions. In both the polarized and unpolarized cases the general formula is 



d(A)a ^ , 1 /•! dv dw ^A, 

dkxidyidz irb 7-^ Jvw I — v Jvw/v w 



)ff{x,,M' 



1 diA)a'\^^^ _ ^^^^^ M\ Ml) 



V dv 



2-n 



(2.17) 
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The indices i,j run over quark flavours or represent a gluon, and in the case of qq scattering, 
the first term is the LO contribution, which is absent for the qg initiated process. 

III. MONTE CARLO CALCULATION 

The combination of analytic and Monte Carlo techniques used here to perform the phase 
space integrals has been documented and described in detail elsewhere (see Bailey et al, 
||T0| and references therein), so our discussion will be fairly brief, highlighting mostly those 
features which are important to our calculation. The basic technique comprises isolating 
those regions of phase space where soft and collinear singularities occur and integrating over 
them analytically in 4 — 2e dimensions. In this way the singularities are exposed as poles 
in e. These regions are isolated form the rest of the 3-body phase space by the imposition 
of arbitrary boundaries between them, achieved by introducing cut-off parameters. The 
soft gluon region of phase space is defined to be the region where the gluon energy, in a 
specified reference frame, usually the subprocess rest frame, falls below a certain threshold, 
5s\/^/2, where 6s is the arbitrary cut-off parameter and s is the center-of-mass energy in 
the parton-parton system. If we label the momenta for the general 3-body process by 
Pi + P2 ^ Ps + Pi + P5, then we can define the general invariants by Sij = {pi + PjY and 
tij = {pi — pjY. The collinear region is then defined as the region where the value of an 
invariant falls below the value 5cS. 

The phase space integration over the mutually exclusive soft and collinear regions are 
performed not on the full matrix elements but on approximate versions of them defined in 
specific ways. In the soft gluon case they are obtained by setting the gluon energy to zero 
everywhere where it occurs in the matrix elements, except in the denominators. This is the 
soft gluon approximation. Similarly, in the collinear case, each invariant which vanishes is 
in turn set to zero everywhere except in the denominator. This is the leading pole approxi- 
mation. The phase space integrals are then performed on these and only the logarithms of 
the cut-off parameters are retained. All positive powers of the cut-off parameters are set to 
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zero. The meaning of all this is that for the method to work the parameters must be kept 
small, otherwise the approximations made will no longer be valid and the method would 
fail. 

Once the phase space integrals have been performed and the soft and coUinear poles are 
exposed, then the virtual contributions, if any, are added to them, at which time all double 
poles and single poles of soft (IR) origin automatically cancel. The remaining coUinear poles 
are then factorized in the parton distribution and fragmentation functions in the usual way, 
at some scale and using some specific factorization scheme. In our case we use the MS 
scheme. At this point one is left with a set of 2-body processes which depend explicitly on 
In 5s and In 5c, and a set of 3-body matrix elements which when integrated over the phase 
space using Monte Carlo techniques, have an implicit dependence on these same logarithms, 
but which now have opposite signs in order to cancel the dependence in the 2-body part. 
The physical cross sections will then be independent of these arbitrary cut-off parameters. 

It is generally very simple to impose experimental cuts and to calculate different observ- 
ables when one is dealing with 2-body matrix elements only, as in the usual LO calculations. 
This is not the case when we consider 3-body processes, as there is the need now to cancel 
soft and coUinear divergences, but beyond that, the standard techniques require often com- 
plex Jacobian transformations to calculate cross sections differential in different variables, 
and the phase space integrals can sometimes only be done analytically when specific lim- 
its of integration are involved. This can sometimes make in impossible to impose cuts on 
the kinematic variables. Thus, fully analytic methods of performing calculations for physi- 
cal processes, although in some cases desirable, can be rather restricted in their usefulness 
when confronted with experimental situations where it is often desirable and sometimes even 
unavoidable that cuts be made. The combined analytic and Monte Carlo method does not 
suffer from these particular drawbacks, in that it is very easy to calculate cross sections 
differential in many different variables at once, and cuts on the kinematic variables to match 
those made in the experiments can be easily imposed because the phase space integrals are 
performed numerically after all singularities have been dealt with. 
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Our calculation of double prompt photon production proceeds along exactly the same 
lines to that described in the second of Ref. |]10[. In fact the only difference is that in the 
polarized case we must deal with 'hat' momentum integrals as outlined in section II. It 
turns out that the 'hat' momenta only contribute in the collinear limit when an initial state 
parton splits into two partons, one of which participates in the hard scattering. This point 
is discussed in fair detail in Ref. [§, where it was suggested that, since these contributions 
seem to be of the same form for a particular vertex in the collinear limit, they may be 
regarded as universal properties of these parton legs in the HVBM scheme, and as such may 



be factorized into the parton distributions under a new (MSp) factorization scheme. In our 
calculation, we choose not to use this scheme because the new polarized parton distributions 



are being calculated in the standard MS scheme. Thus, these contributions are included in 
our 2-body part on the cross section. 



Thus all the discussion in Ref. [|T0[ on the Monte Carlo method can be carried directly 
over to the polarized case, expect that once factorization has been performed the expression 
for the remnants of the hard collinear singularities is now given for the polarized case by 



-{qq 77) 



dv 



2n dv 



X 



+ Af^{xuM') -f^{x,/z,M')AP,,{z) 

Jx2 Z 
Jxi Z 



■-f^{x,/z,M')AP,,{z) 



+ A/f(x2,M2) / -ff{x,/z,M')AP,,{z) 



x-\ Z 



(3.1) 



with 



AP,,(^) = AP,,(z) In (^^c^) - AP:^. 



The polarized Altarelli-Parisi splitting functions in 4 — 2e dimensions are 

'1 + ^2 



(3.2) 



APq,{z,e) = CF 



1 - z 



+ 3e 1 - z 



(3.3) 
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and 

AP,,(z,6) = ^[(2^-l)-2e(l-^)]. (3.4) 

Part of the second terms in the sphtting functions come from the 'hat '-momenta as discussed 
in Ref. 0. The functions AP/^ are defined by the relation, 

AP,,(^, e) = AP,^{z) + eAP;Az). (3.5) 



IV. NUMERICAL RESULTS 

In this section we present some numerical results for polarized and unpolarized double 
prompt photon production at RHIC center-of-mass energies. We first study the non-isolated 
cross section, since we do not know what isolation restrictions will be necessary at RHIC. 
We will then present some results for the isolated cross section, assuming some plausible 
isolation parameters, based on those used by the CDF collaboration at Fermilab. 



Throughout we use the GRV |[T^ parton distributions for the proton in the unpolarized 
case. We use Ag^^ = 0.200, to match the GRV parton distributions, but we do not include 
any contribution for charm quarks, since these are not included in the polarized parton 
distributions. For the polarized parton distributions we use the two sets proposed by Cheng 
and Wai In the first set which we will refer to as scenario a, a large polarized gluon 
distribution is assumed, and the SU(3) flavour symmetric sea quark distribution is assumed 
to vanish at the input scale. In the other case (scenario b), the gluon distribution vanishes at 
the input scale, but the SU(3) flavour symmetric polarized sea quark distribution is assumed 
to be directly related to the unpolarized strange sea, leading to a large negatively polarized 
sea, and small polarized gluon distributions. In both cases the valence distributions are 
assumed proportional to the unpolarized ones. The main aim of this paper is not to provide 
the most up to date phenomenological study of the polarized cross section, but to gauge the 
importance of higher order corrections and make contact with the results of Ref. [12]. We are 
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also interested in whether the cross section is sensitive to AG, hence the use of the Cheng 
and Wai distributions which have extreme gluon distributions. 

For the electromagnetic coupling constant we use aem = 1/137 and we use the two 
loop expression for Q;s(yU^). Unless otherwise stated, we set all factorization/renormalization 
scales to /i^ = {{kj^iY + {kj.2Y)l'^- We make use of the asymptotic parametrizations for the 
parton to photon fragmentation functions provided in Ref. [Q. We assume the maximum 
= 500 GeV for the RHIC centre-of-mass energy. 



In a future study |]T2| we will use photon fragmentation functions evolved in NLO, po- 
larized parton distributions evolved in NLO as well as including a contribution from charm 
quarks. There we hope to make up-to-date estimates of the cross section ate RHIC. 

A. The Inclusive Cross Section 

In the analytic calculation of the inclusive cross section, we exposed the various soft and 
collinear poles by expanding the integrated matrix elements as plus-distributions in the vari- 
able z. The details are given in the Appendix. This procedure ensures that these integrable 
singularities can be treated numerically, but it prevents us from providing distributions in 
the 2;-variable with infinitely sharp resolution. That is, in order for us to present finite z- 
distributions in we must integrate over some finite range of z. Following the procedure used 



in Ref. [1^, we define 

d{A)a 1 d{A)a 



-dz', (4.1) 



dkj'idyidz Az Jz-^ dhj^^dyidz' 

and provide the distributions in z with finite bin widths Az. For the /ct distributions we 
integrate over a specified range of z, 

= n Ji^dz. (4.2) 
dkj^^dyi Jza dk'^^dyidz 

In fig. 4a we show the z-distribution for the unpolarized cross section, at fc^^ = 5 GeV 
and yi = 0, for a bin size Az = 0.2. The LO Born term gives a sharp peak at z = 1 but the 
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effect of the higher order corrections is to reduce and broaden this peak. The higher order 
correction terms, neglecting the gg process which is positive and contributes only at z = 1, 
give a negative contribution at 2; = 1, as indicated by the dashed line. It turns out that 
the fragmentation terms where the trigger photon is produced via fragmentation give large 
contributions in the region -2 < 1, whereas those where the other photon is produced in this 
way give a positive but not as large contribution aX z > 1. This is partly due to the fact 
that phase space runs out in this region. The overall shape of the distribution is similar to 



that obtained in Ref. 10 



In Fig. 4b we compare the polarized and unpolarized cross sections with the same pa- 
rameters as in fig. 4a. The unpolarized cross section is scaled by a factor of 1/10 for easier 
comparison. In general the curves have a similar shape, but due to cancellation between 
various parts of the polarized cross section, it has a much more sharply pronounced peak at 
z = l. 

The distribution is displayed in fig. 5a for 0.2 < z < 2.0 and \yi\ < 3. The solid 
curve is the unpolarized cross section, while the dashed and dotted curves are the polarized 
cross sections for scenario a and h respectively. Clearly, given a reasonably large luminosity, 
the cross sections are large enough to measure out to a kx of about 50 GeV say. Also the 
two scenarios give significantly different predictions indicating sensitivity to the polarized 
distributions. The main question is the sensitivity of the cross section to polarization. To 
determine this we plot the longitudinal asymmetry defined by 

cZAo- 



All = (4.3) 

The dotted curve shows the asymmetry predicted for the Born cross section using the parton 
distributions of scenario a. At the hard parton-parton scattering level this asymmetry is 
All = —1, and is thus modified significantly by folding with the parton distributions. Note 
that this is not the asymmetry for the LO cross section, since we have not included the 
contributions from photon fragmentation. The solid and dashed curves give the asymmetry 
of the full higher cross section for scenario a and b respectively. In the case of scenario a, the 
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asymmetry is already about 10% at kr — 15 GeV, and rises to nearly 20% at kr — 50 GeV. 
For scenario b the asymmetry is more modest, varying between 4 and 9% for a similar range 
in kr- Thus, the inclusive cross section is sensitive to the polarized parton distributions 
used, and the asymmetry not too small in some accessible kinematic regions. 

Of course as we have stated before the inclusive cross section we have calculated is 
probably not going to be of much practical use at the RHIC collider due to the need to 
isolate in order to identify the photon signal. Its main usefulness is as a check on the more 
versatile calculation using Monte Carlo techniques which we shall discuss in the next section. 

B. The Isolated Cross Section 

In this section we study the isolated cross section for double prompt photon production 
at RHIC. We keep all distributions and parameters the same as for the inclusive case, but 
in addition a cone size of R = 0.7 and an energy resolution parameter e = 2 GeV/kj^ is 
used. The cone of radius R, with the photon at the centre, is defined in the pseudorapidity- 
azimuthal angle plane (y — 0-plane) by the relation 

R = ^(Ay)2 + {A4>f. 

If hadronic energy greater than a fraction e of the photon energy is observed in the cone 
then the event is rejected. This serves to define the isolated cross section. 

In fig. 6a we display the kx distribution of the full isolated polarized and unpolarized 
cross sections. Note that this curve cannot be directly compared with that in fig. 5a since 
here the cuts are different. In fig. 6a both photons are still allowed to have rapidities in the 
range \y'^\ < 3, but for a fixed /cy of one of the photons which we plot vs cross section, the 
other is allowed to have /ct > 5 GeV. Our predictions indicate a measurable cross section 
out to kr around 30 to 40 GeV, given enough luminosity. Again, in the polarized case, 
sensitivity to the polarized parton distribution chosen is evident from the figure. 

Fig. 6b shows the asymmetries of the cross sections in fig. 6a for scenarios a and b. There 
is a very pronounced difference in the asymmetries for scenarios a vs b, indicating a corre- 
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spending sensitivity to the polarized parton distributions. The implication of these results 
is that discrimination between the two extreme cases presented here should be possible from 
measurement of this cross section with moderately good statistics at RHIC. In fig.6c we 
break down the asymmetry in the various initial state contributing processes for the case 
of scenario a. The gg initiated process, as expected, only gives a significant contribution at 
low kr due to the steeply falling gluon distribution. But we find that for this scenario the 
asymmetry is dominated by the qg initiated process. We did a similar analysis for scenario 
b and found that the contribution from the qq initiated process was largest and that from 
the gg initiated process was completely negligible. The cross section is thus very sensitive 
to the size of AC 

In fig. 7 we show the scale dependence of the isolated cross section, with the same cuts as 
above, for the unpolarized case. We chose three different scales /i^ = n^Ukj^i)^ + (^t2)^)/2, 
where n — 1,1/2 and 2. All factorization/renormalization scales are varied simultaneously. 
There is obviously some dependence of the cross section on the scales chosen, leading to an 
estimated 20% uncertainty in our predictions at kx = 15 GeV. 

Lastly in fig. 8 we show the kx distribution of the cross section for unpolarized isolated 
double prompt photon production for the same cuts chosen above but at a cms energy of 
\/S = 14 TeV appropriate for the LHC collider. The cross section is as expected much 
larger than at lower energies and even at kr — 100 GeV we can expect a cross section of 
more than 0.1 pb. 

V. CONCLUSIONS 

In this paper we calculated the cross section for polarized and unpolarized double prompt 
photon production in pp collisions at cms energies appropriate for the RHIC collider. We 
examined whether the isolated cross section is sensitive to the spin dependent gluon distri- 
butions of the proton, by comparing two extreme (LO) parametrizations of AG. Our results 
indicate that the cross section is large enough to measure, and that it should indeed be 
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sensitive to the polarized gluon distributions leading to the possibility that it may be useful 
to help discriminate between some extreme scenarios for the polarized distributions. Our 
results differ from those reached in a previous ,mostly LO, study in that we find a significant 
contribution to the asymmetry from the qg initiated process, leading to a corresponding 
sensitivity of the cross section to the polarized gluon distribution, AG. This cross section 
should therefore be useful as a supplement for information on the polarized distributions 
gathered from more sensitive sources such as jet or prompt photon production. The main 
drawback of our study is the use of rather outdated polarized parton distributions evolved in 
LO only. We therefore do not claim to have provided NLO estimates for the polarized cross 
section at RHIC, but only an indication as to the size of the cross section and sensitivity to 
AG. 

Indications are that the cross section for double prompt photon production will, as 
expected, be sizable at the LHC, and will thus provide a significant background to Higgs 
searches. 

VI. ACKNOWLEDGEMENTS 

C.C. thanks the Theory Group at Argonne and in particular Alan White for their kind 
hospitality during the final stage of this work. We thank E. L. Berger, G. Bodwin, R. 
Field, A. White, W. Vogelsang and G. Ramsey for comments and discussions on the matter 
presented, and W. Vogelsang for reading the manuscript. We thank P. Ramond for illu- 
minating discussions concerning the role of the anomalies. Finally we express our warm 
gratitude to Derek and Liz Kruk, Deb Petrie and Betsy Herman for encouragement. This 
work supported in part by the U.S. Department of Energy, Division of High Energy Physics, 
Contract W-31-109-ENG-38 and DEFG05-86-ER-40272. 



19 



APPENDIX A: REGULARIZATION SCHEMES 



We use the 75 scheme of Breitenlohner and Maison , which is free of internal incon- 
sistencies in the definition of 75 in D— dimensions. An anticommuting 75 is not compatible 
with Dimensional Regularization |n,T3,?T|. In the t'Hooft-Veltman scheme, systematized 



by Breitenlohner and Maison, an n-dimensional 7^ matrix is split into its 4-dimensional com- 
ponent 7^ and a remaining component 7. Thus 7/^ = 7^ + 7/x. A suitable representation 
of 75 is 

75 = \e^p,8l"l^l^l'. (Al) 

If it is postulated that 75 anticommute with the other 4-dimensional Dirac matrices and 
anticommute with the remaining ones, then all the corresponding algebraic relations can be 
shown to be consistent with dimensional regularization Therefore we define 

757;. + 7;.75 = 
75% - 7/.75 = 0. 

(A2) 

A discussion of the modifications which appear in the integration over the phase space of 
the 3 final states, due to this ansatz, is presented in the next section. 

APPENDIX B: 3-PARTICLE PHASE SPACE FOR POLARIZED SCATTERINGS 



As we have discussed before, the use of the t'Hooft-Veltman regularization [|l4[ intro- 
duces a dependence of the matrix elements on the hat-momenta which requires, in part, a 
modification of the phase space integral which appear in the unpolarized case. 

In the case of unpolarized scattering, the singularities are generated by poles in the 
matrix elements which have the form l/ta, l/us, l/{t^U'i) and similar ones, in multiple 
combinations of them. Multiple poles can be reduced to sums of combinations of double 
poles by using simple identities among all the invariants and by the repeated use of partial 
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fractioning. This is by now a well established procedure. In our case we encounter new terms 
of the form l/t\ and l/u\ and new matrix elements containing typical factors of the form 
^3, ^2, and ^2 ■ ^3 at the numerator. Let's discuss for a moment these last terms containing 
hat-momenta. It is obvious that by a suitable choice of the parametrizations given by the 
sets 1, 2 ,3 and 4, (defined in the next section) we are able to reduce to the ordinary phase 
space result given by ( |B42D all the matrix elements containing scalar products of the form 
Pi ■ kj, ki ■ kj with i = 1, 2, j = 1, 2, 3. Therefore it is possible to set to zero, after taking the 
traces, all the products contain such combinations of hat-momenta. Then, the only matrix 
elements of hat-momenta which are left and which are not set to zero are those containing 
products of the form ka ■ A;;,, with a, b = 1,2. 

Let's consider the 3-particle phase space integral when hat-momenta are present 

PSs = J rk,rk3d''k2S{kl)S{kl)6{kl)6''{pi +p2-k3-k2- ki)h (Bl) 

and let's lump together the momenta k^ and /c2 as follows 

PS3 = I rk^d''ksd''k2d''k235ikl)5ikl)5ikl)5''ipi+p2 -k2- ki) 

S{P1 +P2-h- k2s)6'\ks + k2- fc23)§'. (B2) 



Eq. ( P2|) , once integrated over /C23, gives ( [Bl|) . We single out the invariant mass of the 
pair (1,2) by the relation 

1 = j dm'^S{m^ - kl^) (B3) 

which we insert in (IB 21) to get 



d''kid''k3d''k2d''k23dm'^5{kl)d{kl)S{kl)S''{pi + P2 ~ k^ - ^23) 
5-(fc3 + k2- k23)5im^ - kl,)k2. (B4) 

Notice that by this trick we can factorize a 2-particle phase space 

p^3 = J d''kidm^d''k236{kl)6{kl^ - m2)5"(pi + p^ - k^ - ^23) x PS2 (B5) 
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PS2 = J (rk2S{kl)S{{k23 - k2f)%. (B6) 
We have set k2 = (/c2, ^2), with 

^2 — (1, cos ^3 sin ^2 sin ^1, cos ^2 sin ^i,cos ^1) (B7) 
being the 4-dimensional part of k2- We easily get 



2 

^2 = ^ sin^ 03 sin^ 9^ sin^ 62. (B8) 



Therefore the usual angular integration measure 

n-l-2 

dft^''-^^^ n ^sm^'-^-^eidOi (B9) 
1=1 

is effectively modified to 

3 ra-2 
^ -Q gi^n-/ ^^^^^ ^ -Q sinn-i-2 g^^g^ (g^Q^ 

Z=l i=4 

This integral is evaluated in a special frame. Assuming that /cfg > 0, we sit in the center 
of mass frame of the (1, 2) pair, in which k23 = {m, 0) and k2 = {E2, 1^2) to get 

PS2 = J dl'k25{kl)5{m^ - 2mE2) 

= k^-'in (Bll) 



where 



m = onrTl' 01 F ^1 F sin-2 82 r d93 sin"-=^ 0s (B12) 

z"i [n/z — zj JO Jo Jo 

where 9i and $2 are the only relevant angles which appear in the matrix elements and 
therefore are not integrated. We have displayed also the 63 integral since it is different from 
the unpolarized case. 

After integration over and k23 and a simple covariantization we get 

PSs = J d^k,5{kl)[(p^ +P2- k^)Y'~' mi (B13) 
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The coUinear or infrared poles are therefore isolated in the cm. system of the pair. The 
remaining part of the integral is evaluated in the cm. system of the two incoming partons 
Pi,P2- At this point let's consider the sub-integral 



J El'-^dEssin^'-^eideiis + t + (B14) 



n/2-l 



T[n/2 - 1] 
We have chosen the parameterization 

/ci = i?3(l, cos6'2 sin 6*1, COS 6*1) (B15) 

where the dots indicate n — 3 components which we integrate over in a trivial way. 

We introduce the change of variables {cosdi, E^) —>■ {v,w) with v and w defined as in 
( p. 31) . In the cm system of the two incoming (massless) partons we have that pi = Qn^ 
and p2 = Qn~ with Q = \J s/2 

n± = -^(l,0^,±l). (B16) 

From t = {ki — piY = —^^£'3(1 — cos 6*1), u = {ki—p2Y = —y/sE^^l + cos 6*1) and using 
3D we get 

vw — 1 + V 



COS^^l 



vw + 1 — V 



E3 = ^{vw + l-v) 
s + t + u = sv{l — w) 

(l-cos^^0 = 7^^^^^ (BIT) 

and the jacobian of the transformation to be d{cosOi, E^) / d{v,w) = y/sv/{l — v + vw). 
Therefore we obtain 

n-5/2 ^ 
^ 2"+ir[n/2-l/2]r[n/2-2] J v ; v ; 

X Trf^i rd^2sin"-i^isin"-2^2. (B18) 
Jo Jo 
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In order to compare this result with pO[ we need an additional normalization factor 
l/(27r)^~''% due to the different definition of PS3 we have adopted (now n = 4 — 2e) and we 



get 



32-2.^1-2. / _^ 



24(27r)5-4^ Ul-e) 



2^-2^ / dvdw{l - vy'{l - wf-'v^-^'w-' 



de2sin'-''e2 dOism-'-^'ei. (B19) 
^0 

In order to integrate over the matrix elements, we need to evaluate the various scalar 
products which appear in such matrix elements, in the cm. frame of the pair (1,2). For 
this purpose we define the functions 



P[x,y,z] 



+ y'^ + z'^ — 2xy — 2yz — 2xz^^ 
Ax 



/2 



E[x,y,z] = ^±y^. (B20) 
2vx 



It is easy to show that 



bll = P[s23,Pl,Ui\ 
\P2\ = P[s23,pl,ti] 
Ihl = 1^2! = P[s23,PvPI] 

\ki\ = J — P[s,kl,S23] 
V S23 

Pi = E[s23,pI,Ui] 

pl = E[s23,plh] (B21) 



When all the external lines of the 2 — > 3 process are massless, then we specialize ( P21|) as 
follows 



p", = E[s23,0,Us] 



2^ 

p^ = E[s23,0,t3] 



r-.r_ n J. T ^23 — ^3 



(B22) 



24 



Now, using ts = s{v — 1) and U3 = —svw we get 



.0 _ 



In the massless case, the energies of the external final state particles are given by 



V ■S23 

1 (S - S23) 



^ 2 
1 



s(l - V + w) 

In the derivation of ( P24|) we have used the relation s + ti + ui = S23 together with 



kl = kl = (B24) 



( p.3|) . There are four different parametrizations of the integration momenta which we will 
be using. In the first one, which is suitable for unpolarized scattering [0, one defines (in 
the cm. frame of the pair (1, 2)) 



set 1 



ki = s(l — V + w){l, 0, sin ip, cos ip) 

2vS23 



^3 = 2^ (1; cos^2 sin 9i, COS^i) 



^2 = (1; ~ cos6'2 sin^^i, — cos6'i) 
sv 

Pi = „ , — (1,0, ...,0, sin V^s, cost's) 
2V^ 

P2= \ ^ ^ (l,0,...0,sin^2,cosz^2) (B25) 

2V'523 



where the dots denote the remaining n — 2 polar components. Similarly, in the evaluation 
of the integrals over the hat-momenta we need the other parametrizations 

• set 2 
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Pi =p?(l,0,..., 0,0,1) 

P2 =p2(l,0,...,-sint/'2,0,cost/'2) 
— 0, — sin ■00, 0, cos'^o) 

(B26) 

• set 3 

Pi =p?(l,0, ...,sinV'2,0,cosV'2) 

P2=P^(1,0,..., 0,0,1) 

ki — A;?(l, 0, siii'^i, 0, cos'0i) 

(B27) 

• set 4 

pi =p?(l,0, ...,sinV'o,0,cosV'o) 
P2 =P2((l5 0,...,-sint/'i,0,cost/'i) 
A;i = A;?(l,0,.., 0,0,1) 

(B28) 



where 0, ... refers to n — 5 components identically zero. It is straightforward to obtain the 
relations 



. 2^w{l-v){l-w) 

SUllpQ — - 



1 — w + vw 



2vJw{\-v){l-w) 

^^^^1 = T\ ^7i 7 V 

(1 — vw)\ \ — V + vw) 

sinV^2 = -■ (B29) 

\ — vw 

Which set of parametrizations we are going to use depends on the form of the hat-momenta 
which appear at the numerators of the matrix elements after the traces are performed. 
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Notice that 

k3 + k2 = ^{1,0). (B30) 



In the cm. of the pair (1,2) pi, p2 and ki lie on a plane, with the spatial components 
satisfying the condition Pi + P2 + ki = 0. Using ( P25|) and the expressions of and 
■S23 = sv{l — w) in terms of v and w, we can easily obtain the relations 



sin ip2 



1 — w 



sin ipi = — sin ip 



1 — vw 

'1 — V — vw 



1 — V + vw 



cos^ = . r^-"^ (B31) 
V 1 — vw 



Following Ref. |]TOl we introduce the variable 



where the perpendicular components are measured in the cm. system of the two incoming 
partons pi and p2- We can covariantize ([B32| ) in a trivial manner by using light cone identities 

ki = '^k^k^ — {k^Y = 

{kif = 2ktkz = 2h:Plh:Pl (B33) 

P1-P2 

Similarly, expanding k^ ■ k2 in its light cone components and covariantazing we get 

. k^ = -k, . k2 + k3-P2k.-p, + k,.p,k2-p2 

Pi ■P2 

and 

sks ■ k2 + uk2 ■ pi + tk2 ■ P2 

z = 

tu 

= m-k2, (B35) 



where the four- vector m is defined by 



_ sks + tp2 + upi 

m = . B36 

tu ^ ' 
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( P35|) is the covariant expression of z. At this point, however, it is necessary to evaluate m 
in the cm. system of the pair (1,2). 

Using ( [B25[ ) in Eq. ( [B36| ) it is a simple exercise to show that m has only longitudinal 



components given by 



m = {-y/'{J^ ^,0,...,0, 

tu ^ 1 — w 



1 — vw, 
1-w ' 



= (^)i/2K,0,...,0,m;). (B37) 

Notice that cosip = rriQ/m'^ = tanh x- 

In order to isolate a pair (1,2) of a given z, we introduce the identity 

dz5{z-m.k2) (B38) 
in the expression of the phase space Ps^ which becomes 

2-2e l-2e / _p \ 

^ 24(27r)5-4^28r[l -2e] V(l -e); 7 ^ ' ^ ' 

X de2dz6{z -m-k2) sin^-^' 62 F rf^isin^-^^ O^. (B39) 

The integration over 9i can now be performed and the remaining ^-function eliminated. 
One gets 

r dOi sin^-^^ OMz - m ■ ^2) 
Jo 

= r dOi sin^~^' Oidiz -1/2 + 1/2 cos Oi coth x) 
Jo 

= 2 tanh X9iv, w, z)^"' (B40) 

where 

, , l-w + 4:w{l-v)z{l- z) 

g{v, w, z) = B41 

1 — vw 

The final phase space therefore can be cast in the form 
X r sin^"^^ ^2 

JO 

(B42) 
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APPENDIX C: EVALUATION OF THE PHASE SPACE INTEGRALS 



Prom the definition of the hypergeometric function 

''''^^"r[6]r[c-6]io {l-z/2 + z/2coser ^ ^ 



we easily get 

^l"' '"-'^^L (a^Pcoser = 5=2l^^l«/2- a/2 + 1/2,6+1/2, f^W] (C2) 
We now use the relation 

F[a, P, 7, ^] = (1 - ^)^-"-^F[7 - a, 7 - /3, 7, z] (C3) 

to get 

^2. 2 - 2.1 = gg pp_^',l;'„;!^;,.).,.,. f[l - . 1/2 -e,2-e, (C4) 

and 

At the end the result is expressed in terms of "plus" distributions using various identities 
whose derivation is briefly discussed in the next appendix. 

For instance, let's denote by Ps^lts] the phase space contribution due to a factor l/t^ in 
the matrix element. We get 

Pssih] = ^K^n, p— ^^F[l/2 - e, -e, 1 - e, (3^/a\ (C6) 

where the factor K^np is the same as calculated in the unpolarized case 

Notice that at 2; = 1, we have (5 — a and Taylor expanding the hypergeometric function 
around 2; = 1 we get (with r{z) = P'^/a^) 
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a h 



F[a, 6, c, r{z)] = F[a, b, c, 1] + — ^[a + 1, & + 1, c + 1, l]r'(l 
= F[a,6,c, 1] + (z- l)O(e). 



(C8) 



Then we can set F[a,b,c,r{z)] = F[a,b,c,l] + 0(e) with F[a,b,c,l] = 1/2^' + 0{e^). We 
finally get (see appendix C) 



'1-z) eiz-i 



1 



unp 



'l-z) 



+ r — TT^ + '^(^-i)( iog^> 



(C9) 



''30 VV-^ ^ + 

As an example of applications of of the methods discussed above in the polarized case, let's 
consider the contribution to the final phase space coming from matrix elements of the form 
k /tl We get 



k 



^''^ tl I ~tio^"'"V2(l-e)/ \9 J \l-z\^+^' 



1 



where 



■r]{z)g{z)F[l-e,l/2-e,2-e,(3ya^' 



(CIO) 



viz) 



(l + tanh^x(l-22)) 



(Cll) 



and 



. = — - — — -^v^~^'w~'(l - wy-'il - vY 



(C12) 



For future purposes it is convenient to introduce the function <t{z) ^ q{z)g{z). We get 



2 -"^Po' 



"-30 



1 - e 



<5(z-l) 



(C13) 



M30 VI -e/ 



(C14) 



Notice that in the unpolarized case - as discussed in Ref. - in some specific matrix 
elements, such as PS3 (1/(^13^3)), singularities in both variables v and w are encountered 
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(for w = 1 and z = 1). The regularization in terms of plus-functions of the corresponding 
contributions and a detailed discussion of the derivation, which is similar to the polarized 
case can be found in Ref. UTOl . 



APPENDIX D: IDENTITIES FOR "PLUS" DISTRIBUTION 



As we have discussed before, we take z to be positive from the beginning. 

In the expression of the final cross section, we are going to encounter singularities at 
z = 1 (and at w = 1) which have to be regulated in an appropriate way by the use of "plus" 
distribution. Identities for "plus" distribution can be easily derived by the integration by 
parts method, as we are going to illustrate briefiy by an example. Assuming 2; > 0, we define 
two kinds of regulated "plus' distributions 9{1 — z)/{l — 2)+ and 9{z — l)/{z — 1)+ by 



(1 — ZjA. Jo 1 — z 



(Dl) 



Let's define the two functions 



h{z) 



1 + tanh^ x(l - 2^) 

2 tanh x 
I - w + Azw){l - v){l - z) 
1 — vw 



(D2) 



with h{z„ 



dzf[z 

+m (- 



Zmin- Then we get the identities 
9iz-l) fh^[z]Y 



[Z - 1)1+2^ 
1 



9[z] 



dzf[z] 



e{z-i) fh^[z] 



(z - iy+'^ \ g[z] 



2e 2 g^[l\ 



(D3) 



We have integrated once by parts and the boundary (regular) terms have been set to zero 
before taking the limit e — > and expanding in e. The procedure has to be repeated once 
more if the boundary terms are still singular at the edge of the interval of integration. Now 
we use Zjnax — 1 = —Zmin together with the simple identity 
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e{z-i) 



g[z] J {Z - 

to finally obtain the relation 



e[z - 1] 



(z-iy 



0[e] 



(D4) 



9{z-l) 



'h^[z]Y 1 
^g[z]) {z-iy-^ 



9[z - 1] 



(^-1). 



+ S{z-1) 



-— + log(-z„ 



1 



(D5) 



In order to derive a similar identity in the case of < 2; < 1 we proceed in a similar way 



dzf[z] 



(z-l) 



l+2e 



9[z] 



dzf[z] 



{Z - 1)1+2. ^ ^[^] 



Further manipulations similar to those presented above then give 



^[1 - z 



1 



g[z] {z-iy+ 



9[1 - z] 



1 



(^-1) 

Combining (^^) and (|D7|) , after some manipulations we get the identity 
1 fh^[z]\ 0(1 -z) eiz-1) ^.^ X, 1 , 



In a similar way we can derive the identity 



'h\zy 



r]iz)giz)F[z] = 9il-z) 



1 - z 



\z - \ g{z) 

+6{z-l)[-^r^{z)g{z)^ -U{z - g{l) F[l]+ 0(1) 



(D6) 



(D7) 



(D8) 



(D9) 



Using F{1] = 2(1 + 0(e)) where F[z] = F[l - e, 1/2 - e, 2 - e, r{z)] we get ( |CT3|) and (|CT|). 
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Figure Captions 

[1] (a) Lowest order Feynman diagrams for double prompt photon production, (b) 
Examples of single fragmentation contributions to double prompt photon pro- 
duction, (c) Example of box diagram contribution to double prompt photon 
production. 

[2] (a) Examples of virtual corrections to the lowest order diagrams, (b) Examples 
of next-to-leading order three-body final-state diagrams for the qq initial state. 

[3] Examples of contributions to the higher order qg initiated process. 

[4] Unpolarized non-isolated cross section da / dkj'idyidz as a function of z for p+p — > 
+ + -'^ at ^/s = 500 GeV. We set yi — 0. Results are presented in the form 
of a histogram in bins of width Az = 0.2. In (a), for kj^^ = 5 GeV, we show 
the net contribution from the lowest order process qq — > 77 and from all the 
higher order processes and the full sum. In (b) we compare the polarized and 
unpolarized cross sections for the same parameters as in (a). The unpolarized 
has been multiplied by 0.1 for easier comparison. 

[5] (a) The transverse momentum dependence of da/dkj'idyidz, the non-isolated 
cross section, for z integrated over the interval 0.2 < 2; < 2.0. The upper solid 
line shows the full unpolarized cross section. The lower curves show the polarized 
cross section as given by the two different scenarios for the polarized parton 
distributions discussed in the text, (b) The longitudinal asymmetry, defined in 
the text, for the non- isolated cross section as predicted by parton distributions 
assuming scenarios a and b. The asymmetry for the lowest order Born cross 
section is included for comparison. 

[6] (a) kr dependence of the isolated cross section integrated over rapidity range 
— 3 < y < 3 using isolation parameters given in the text. The unpolarized cross 
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section as well as the polarized cross section assuming scenarios a and b. (b) 
The longitudinal asymmetry for the cross section in (a) as predicted by scenario 
a and b. (c) The same asymmetry as in (b), but only for scenario a broken down 
into contributions from the qq, qg and gg initiated processes. 

[7] The renormalization/factorization scale // dependence. For the sum of all con- 
tributing subprocesses, da/dkj^^, for — 3 < y < 3 is shown as a function of /cj^ 
for three values of = + (/cj.2)^)/2: 0.5, 1.0, and 2. 

[8] The same cross section as in fig. 6a for the unpolarized case but with y/S — 14 
TeV. 
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